Self force via a Green's function decomposition 
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The gravitational field in a neighborhood of a particle of small mass fi moving through curved 
spacetime is naturally decomposed into two parts each of which satisfies the perturbed Einstein 
equations through O(n). One part is an inhomogeneous field which looks like the n/r field tidally 
distorted by the local Riemann tensor. The other part is a homogeneous field that completely 
determines the self force of the particle interacting with its own gravitational field, which changes 
the worldline at 0{fi) and includes the effects of radiation reaction. Surprisingly, a local observer 
measuring the gravitational field in a neighborhood of a freely moving particle sees geodesic motion 
of the particle in a perturbed vacuum geometry and would be unaware of the existence of radiation 
at O(fi). In the light of all previous work this is quite an unexpected result. 

PACS numbers: 04.25.-g, 04.20.-q, 04.25.Nx, 04.30.Db 



INTRODUCTION 



The principle of equivalence in General Relativity im- 
plies that a freely moving particle, of infinitesimal mass 
and size, follows a geodesic of spacetime [Q. For a small 
but finite mass /i, the particle perturbs the spacetime ge- 
ometry at 0(/i). The interaction of the particle with this 
perturbation of the metric is called the "self force" which 
changes the worldline at 0{n) and includes the effects of 
radiation reaction. 

A local observer measuring the gravitational field in 
a neighborhood of the particle, with no a priori knowl- 
edge of the background geometry, sees a combination of 
the background metric plus its 0{fi) perturbation caused 
by the particle. Local measurements, in a neighbor- 
hood of the particle, cannot distinguish these two spe- 
cific contributions. Nevertheless, we show below that the 
combined metric in the neighborhood of the particle can 
be uniquely decomposed into two distinct parts: (i) the 
tidally distorted n/r piece of the perturbation of the met- 
ric, which is an inhomogeneous solution of the perturbed 
Einstein equations at 0(/i), and {ii) the sum of the back- 
ground metric and the remainder of the metric perturba- 
tion, which together arc a homogeneous solution of the 
Einstein equations through 0(/i). 

The self force is shown to result in geodesic motion 
of the particle, through 0(/i), in the homogeneous per- 
turbed metric, {ii) above. With only local measurements, 
the observer has no means of distinguishing the homoge- 
neous perturbation from the background metric through 
0{n). As the particle moves along a geodesic of this 
perturbed vacuum geometry, the observer sees no effect 
which one would be compelled to interpret as radiation 
reaction at 0(/i). In this precisely defined fashion, we ex- 
tend the understanding of the principle of equivalence in 
General Relativity: The orbit of the particle is a geodesic 
through 0{/i), in a locally measurable, vacuum gravita- 
tional field. 



II. HISTORICAL BACKGROUND 

In Dirac's ||] analysis of the self force for the electro- 
magnetic field of a particle in flat spacetime, he used the 
conservation of the stress-energy tensor inside a narrow 
world tube surrounding the particle's worldline to derive 
equations of motion with radiation reaction effects in- 
cluded. DeWitt and Brehme H extended this approach 
to allow for curvature in spacetime. Mino, Sasaki and 
Tanaka @ generalized it to include the gravitational self 
force as well. Quinn and Wald(5) and Quinn|| have ob- 
tained similar results by pursuing an axiomatic approach 
for the gravitational, electromagnetic and scalar field self 
forces. 

On a formal level, these analyses provide a coherent 
and consistent treatment of the self force. Analysis be- 
gins with a given worldline T, described by z a {r), which 
obeys the Lorentz force law through the fixed background 
gravitational and electromagnetic fields. The self force 
gives the actual worldline an acceleration away from T, 
and it is the mass times this acceleration which it is 
wished to determine. 

Towards this end, first the retarded field is obtained 
in terms of the corresponding Green's function. Tradi- 
tionally the Green's function has been decomposed into 
"direct" and "tail" parts. The resulting "direct" part of 
the field at a point x is determined completely by sources 
on the past null cone of a;; in flat spacetime this is the 
Lienard-Wiechert potential for elcctromagnetism. The 
curvature of spacetime allows for an additional contribu- 
tion from sources within the past null cone, and this is the 
"tail" part. The self force on a particle is then composed 
of two pieces: The first piece comes from the direct part 
of the field and the acceleration of T in the background 
geometry; in flat spacetime this is the Abraham-Lorentz- 
Dirac (ALD) force. The second piece comes from the tail 
part of the field and is present in curved space even if T is 
a geodesic. While the decomposition of the field into the 
direct and tail parts has been useful for describing the 
self force, neither of these parts individually is a solution 
of the field equation. 
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In flat spacetime Dirac decomposed the retarded 
electromagnetic field into two parts: (i) the "mean of 
the advanced and retarded fields" which is an inhomo- 
geneous field resembling the Coulomb q/r piece of the 
scalar potential near the particle, and (ii) a "radiation" 
field, his Eqs. (11) and (13), which is a homogeneous so- 
lution of Maxwell's equations. He described the self force 
as the interaction of the particle with the radiation field, 
a well-defined vacuum Maxwell field. 

Previous descriptions of the self force in curved space- 
time H [|, ||, H reduce to Dirac's result in the flat space- 
time limit. And they provide clear, adequate expressions 
for the self force. However, they do not share the physical 
simplicity of Dirac's analysis where the force is described 
entirely in terms of an identifiable, vacuum solution of the 
field equations. Indeed, for an electromagnetic field, the 
vector potential A s ° u [cf. Eq. ( fL4| ) for the corresponding 
scalar field expression] might be said to be responsible for 
the self force. But generally A® elf does not satisfy Eq. fl23| ) 
below, and the current density J a , which would result 
from th e ap plication of the operator on the left hand side 
of Eq. (H) to A s a Ql{ , would have a non-van shing charge 
distribution in the vicinity of the particle. Furthermore, 
A 11 ® is non-differentiable at the location of the particle 
with four-velocity u a if (Rab — \g a bR)u h ^ there. Sim- 
ilar statements hold for scalar and gravitational fields. 

In this letter we present a curved space generalization 
of Dirac's scheme. We find that, near the worldline F, 
even in a curved geometry, the retarded field can be de- 
composed into two parts such that the first, A^, is an in- 
homogeneous solution to the field equation with a point 
source, similar to the "mean of the advanced and re- 
tarded fields," while the second, A^, is a homogeneous 
solution which yields a complete description of the self 
force. These parts are related to but distinct from the 
usual direct and tail parts. Scalar and gravitational fields 
are analyzed similarly. 

Below, we first focus on the scalar field case, and then 
briefly describe our results for the electromagnetic and 
gravitational fields. 



III. SCALAR FIELD 

The force on a point particle with charge q moving in 
a scalar field may be deduced from (see discussion in M ) 



F a = qV a ip, 



(1) 



where the derivative of the field is to be evaluated at the 
location of the particle. Usually it is implicit that ip does 
not include the field from the particle itself but is com- 
posed only of the "external" field. The description of the 
self force is not so straightforward. However, one can in- 
troduce a quantity i/j seli , as in our Eq. ([l4|), which when 
substituted into the right hand side of Eq. (Q) formally 
yields the scalar-field self force as described by Quinn Q . 
Similar expressions for the force in terms of the deriva- 



tives of the field are given in Eq. (^8|) and Eq. (|34|) for 
electromagnetic and gravitational fields. 
The scalar field equation 

V 2 V> = V V a V> = -4irp (2) 

is formally solved in terms of a Green's function, 

V 2 G(x,z) = -(-gr^ 2 S 4 (x-z). (3) 

The source function for a point charge moving along a 
worldline T, described by coordinates z a (r), is 



p{x) 



(-gy^S^x-z^dT, 



(4) 



where r is the proper time along the worldline of the 
particle with scalar charge q. The scalar field of this 
particle is 



ip(x) = Anq I G[x,z(t)]& 



(5) 



A symmetric scalar field Green's function is derived 
from the Hadamard form to be 

G sym (x, z) = ^- [u(x, z)S(a) - v(x, z)Q(-a)\ , (6) 

where u(x, z) and v{x, z) are bi-scalars, the properties of 
which are described by DeWitt and Brehme They are 
determined by a local expansion in the vicinity of T, and 
are symmetric under interchange of x and z. a is half of 
the square of the distance measured along the geodesic 
from x to z with a < for a timelike geodesic, and a = 
on the past and future null cones of x. The expansions 
for the bi-scalars u(x, z) and v(x, z) are known to be con- 
vergent within a finite neighborhood of T if the geometry 
is analytic J?]]. The 0(— a) guarantees that only when x 
and z are timelike related is there a contribution from 
v{x, z). The terms in any Green's function containing u 
and v are frequently referred to as the "direct" and "tail" 
parts, respectively. 

The direct part of G sym (x, z), has support only on the 
null cone of x, and the resultant direct part of the field 
for a particle moving along T is 



qu(x, z) 


+ 


qu(x, z) 


2(7 


Tret 


2& 



(7) 



where a = da(x, z(r))/dr, and r ro t and r a( jv refer to the 
proper time of the intersection of T with the past and 
future null cones of x, respectively. 

The tail part of G sym (:r,z), has support within both 
the past and future null cones of x, and the resultant tail 
part of the field for a particle moving along T is 



v(x,z)dr. (8) 



DeWitt and Brehme g use local expansions in the 
vicinity of T to show that |)| 

1 



u{x, z) = l + —R ab V a a W + Q{r d ), 



(9) 
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r is the proper distance from x to T measured along the 
spatial geodesic which is orthogonal to T. They also show 
that the symmetric bi-scalar v(x, z) is a solution of the 
homogeneous wave equation, 



V 2 v{x,z) = 0, 



and that for x close to T, 
1 



v(x, z) 
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R{z) + 0(r), x^T. 



(10) 



(11) 



The retarded and advanced Green's functions are de- 
rived from G sym (x, z), 



G" ct {x,z) 
G &Av {x,z) 



2Q[H{x),z]G syai {x,z) 
2Q[ Zl T,(x)]G sym (x,z), (12) 



where <d[S(x), z] = 1 — 8[z, £(x)] equals 1 if z is in the 
past of a spacelike hypersurface S(x) that intersects x, 
and is zero otherwise. The G lad (x,z) implicitly used by 
Dirac in flat spacetime is 



G rad (a;, z) = G rot (x, z) - G sym (x, z). 



(13) 



Note that the fields resulting from two different Green's 
functions that each obey Eq. (|3j) necessarily differ by a 
homogeneous solution of Eq. (0) . 

G Tet (x,z) has reasonable causal structure, and we as- 
sume for simplicity that ip rct is, in fact, the actual field 
resulting from the source particle. 



IV. SELF FORCE 

Careful analyses j| |J, ||, || show that contributions to 
the self force result from both the source's acceleration, if 
r is not a geodesic, and from the curvature of spacetime. 
For these two distinct possibilities, the self force is a con- 
sequence of the particle interacting with either the direct 
or the tail part of its field, respectively. By following the 
detailed derivations in [|| ||, [|, [J, we see that the self 
force may be considered to result, via Eq. ([!]), from the 
interaction of the particle with the quantity 



if, 



self . 



qu(x, z) 


Tadv /•' 


2& 


Tret J ~ 



v[x, z(r)]dr (14) 



where, unlike Dirac's radiation field, ip is not a homo- 
geneous solution of the field equation. 

The first expression in Eq. (^) is finite and differen- 
tiable in the coincidence limit, x — > T. When substituted 
into the right hand side of Eq. (Q) this expression pro- 
vides the curved space generalization of the ALD force, 
and local expansions of u(x, z) and &(x, z) in ^ [|, ||, ||] 
give the resultant force in terms of the acceleration of T 
and components of the Riemann tensor. 

The integral in Eq. ( |l4] ) comes from the tail part of the 
Green's function. Its derivative results, in part, from the 



implicit dependence of r rct upon x. Quinn 
this contribution to the derivative to be 

- qv[x, z(r rot )]V a T rct = q [v&~ l T7 a a\ 

L J ' ret 

qR{x) 



computes 



12r 



(i. - z a ) + 0(r), x^T, (15) 



the spatial part of which is not defined when x is on T. 
In the usual self force analysis, one first averages V a ^ se 
over a small, spatial two-sphere surrounding the parti- 
cle, thus removing the spatial part of Eq. (|l5|). Then 
one takes the limit of this average as the radius of the 
two-sphere goes to zero, thereby obtaining a well defined 
contribution to the self force pi 0, 0, [I. 

V. HOMOGENEOUS FIELD 

We now provide an alternative expression for the field 
responsible for the self force. 

Given one Green's function which solves Eq. (||), a sec- 
ond may be generated by adding to the first any bi-scalar 
which is a homogeneous solution of Eq. (^) . v(x, z) is just 
such a bi-scalar and is also symmetric, v(x, z) — v(z, x) 
Thus, a new symmetric Green's function is 

G s (x,z) = G sym (x,z) + ^-v(x,z) 



8tt 
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= —[u(x,z)5(a)+v(x,z)G{a)], (16) 

which has no support within the null cone. We use 
G s (x, z) only in a local neighborhood of the particle, and 
do not depend upon any knowledge of its global existence. 
G s (x, z) does have support on the null cone of x, just as 
G sym (a;, z) does, and also outside the null cone, at space- 
like separated points. The use of G s (x, z) is thus not 
complicated by the need for knowledge of the entire past 
history of the source and is amenable to local analysis. 
The corresponding field 



i>*{x) 



qu(x, z) 




qu(x, z) 


2& 


+ 

Tret 


2a 



v(x, z) dr 



(17) 



is an inhomogeneous solution of Eq. (Q) just as ip ret is. 
In the pioneering spirit of Dirac, it is natural to define 



G K (x, z) ee G rct (x, z) - G s (x, z) 



(18) 



[cf. Eq. Remarkably, like G rot (x,z), G R (x,z) has 

no support inside the future null cone. Corresponding to 
G R (x,z), we construct 



qu{x, z) 



2a 



Tret 

1 r T ^ 



v{x,z)dr (19) 
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analogous to Dirac's radiation field. By construction 
is a homogeneous solution of Eq. (||) and has the conse- 
quent property that it is smooth in the coincidence limit, 
x — > r. We note its relation to i/> sclf 



self 



1 



v(x, z) dr. 



(20) 



Our most significant technical result is, in fact, that 
the self force is determined by the particle's interaction 
with , R , since the integral term in Eq. ( p0| ) gives no 
contribution to a self force. For a field point x which is 
near T, the integrand may be expanded using Eq. (11). 
The dominant part of the integral from r re t to r ac i v then 
brings in a factor ofr a dv — T ro t = 2r + 0(r 2 ) times v(x,x), 
and the integral term of Eq. ( p(i| ) is 

— qrv(x, x) + 0(r 2 ) = —qrR(x) + 0(r 2 ), x -> T. (21) 

The derivative of j^qrR{x) gives an outward pointing, 
spatial unit vector near T; this exactly cancels the trou- 
blesome part of V a -0 self in Eq. (15) which is thus absent 
from V a ^ R . The derivative of the remainder term 0(r 2 ) 
is zero in the limit that x approaches Y and gives no 
contribution to the self force. Thus, the self force may 
be seen to be due exclusively to the interaction of the 
particle with ?/> R via Eq. (|l]). We regard this approach 
as preferable, because ip R is diffcrcntiablc at the location 
of the particle, so that averaging is no longer required in 
computing the self force. Even more importantly, ?/? R is 
a homogeneous solution of Eq. . 



VI. ELECTROMAGNETIC AND 
GRAVITATIONAL FIELDS 

The analysis for the scalar field is easily generalizable 
to both electromagnetic and gravitational fields by the 
addition of extra indices to u(x, z) and v[x, z) to create 
corresponding bi-vectors and bi-tensors, and by the in- 
troduction of g ab (x, z), which is the bi- vector of parallel 
displacement ||. The primed indices below refer to the 
source point z, the unprimed indices to the field point x 
as above. The definitions and relationships for the var- 
ious Green's functions follow the same pattern as above 
and are not repeated below. 

For the electromagnetic field, the Lorentz gauge re- 
quires 



VaA U = 0, 

so that Maxwell's equations become 

v 2 A a _ R a bA b = _ 47r ja < 

DeWitt and Brehme H show that 

u ab , (x, z) = g aV (x, z)u(x, z), 



(22) 



-^9b'c'R 



)+0(r). (25) 



and that (as x — > T) 

v ab >(x,z) = \^g a c {r W{ 

If (R a b — \gabR)u b ^ at the particle, then ^4„ elf is non- 
differentiable there. 

Similar to the scalar field, is an inhomogeneous 
solution of Eq. (|2^) . That it also satisfies the gauge con- 
dition ( |22| ) follows from an argument similar to that after 
Eqn. (3.37) in g. We have 



V a ^ = Jv a ,(G s J a ')V=9d 4 x', 



(26) 



assuming that J a is conserved, where G s is the scalar 
Green's function of Eq. (fut). That G s has no support 
within the past or future null cone implies that the inte- 
gral is zero when written as a boundary integral. 
With the definition 



A, 



Al 



A„ 



(27) 

A R gives a homogeneous electromagnetic field. The elec- 
tromagnetic self force becomes 



F a = qg ac (V c Af - V b Af)z b . 



(28) 



This combines with the Lorentz force law from the back- 
ground to determine the actual worldline of the particle. 
For the gravitational field, the harmonic gauge requires 

V a h ab = 0, (29) 

and, with R ab — 0, the perturbed Einstein equations are 

h a b — kgabh c r . is the trace reversed version 



2R a C b d h cd = -167rT ab , (30) 



where h ab = h ab - |i 
of the metric perturbation h ab . Mino, Sasaki and Tanaka 
fl show that 



u a bc'd'{x,z) = 2g ac ,(x,z)g bd >{x,z)u(x,z), (31) 

and that (as x — > T) 

v a bc'd>{x, z) = -g a e g b f Rc'e'd'f(z) + 0(r). (32) 

If R ca dbU c u d ^ at the particle, then h s * h is non- 
differentiable there. 

Similar to the electromagnetic field case, h^ b is an in- 



homogeneous solution of Eq. (30) satisfying the gauge 
condition (E9|), if T ab is conserved. With the definition 



a ab — 11 ab n ab^ 



/i R b is differentiable on T, and 



pa 



fj,(g ab + z a z b )z c z d {V c h 



db 



(33) 



:Vbt&). (34) 



Subject to this force, the particle moves along a worldline 

(23) which is actually a geodesic for a metric composed of 
the background geometry complemented by h^ b . While 
geodesic motion has been demonstrated in the past 

g, only in our case is the reference metric, g ab + h^ b , a 

(24) vacuum solution of the Einstein equations through O(fi). 
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VII. DISCUSSION 

For the clearest demonstration of the impact of our 
analysis, we consider the free motion of a particle of small 
mass n in the purely gravitational case. With no a pri- 
ori knowledge of the background geometry, an observer 
makes local measurements of the metric within a neigh- 
borhood of the worldline of the particle. That field has 
two distinct contributions. The first is the background 
metric combined with h^ b — this combination appears 
as an "external", homogeneous field; no local measure- 
ment distinguishes h^ b from the background. The sec- 
ond comes from ft,^ b ; for free motion the observer knows 
this contribution to be simply the fi/r part of the metric 
with its tidal distortion from the "external," homoge- 
neous field As a consequence of Eq. the local 



observer naturally sees that the worldline of the particle 
is a geodesic in the combined, "external" homogeneous 
field which he measures. Making only local measure- 
ments near the worldline, the observer sees no radiation, 
no local source for the "external" field and also no effect 
which he would be compelled to describe as radiation 
reaction. 
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